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摘 要 致 动 面 模型 利用 无 厚度 的 平面 代 殖 风力 机 叶片 , 在 平面 上 施加 不 连续 压力 来 模拟 叶片 对 气流 的 作用 , 结合 N-S 
方程 , 在 FLUENT 中 进行 数值 模拟 计算 , 是 对 致 动 线 方法 的 延伸 与 改进 . 运用 该 方法 能 够 简化 风力 机 的 模型 ， 从 而 减少 网 
格 数量 和 计算 时 间 . 采用 线性 分 布下 的 致 动 面 模型 ， 一 种 致 动 面 网 格 的 辨识 方法 ,对 单 台 Nibe A 型 风力 机 的 远近 尾 
流 区 域 进行 模拟 计算 , 包括 尾 流风 速 变 化 ， 泣 流 强度 ， 涡 结构 ,并 将 数值 模拟 的 结果 与 致 动 线 模 型 计算 结果 以 及 实验 数据 
进行 对 比分 析 ， EREEREER RHET, pae aiar p ei 性 以 及 用 于 风力 机 尾 流 场 计算 的 可 行 性 . 
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Numerical Study of Wind Turbine Wake Modeling Based on a 
Actuator Surface Model 
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Abstract In the Actuator Surface Model (ALM), the turbine blades are represented by porous 
surfaces of velocity and pressure discontinuities to model the action of lifting surfaces on the flow. 
'The numerical simulation is implemented on FLUENT platform combined with N-S equations. This 
model is improved on the basis of actuator line model(ALM). By using ASM, the model of turbine 
can be simplified and the quantity of grids and computing time can be significantly reduced. A linear 
distribution model and a ASM Grid identification method are presented. This paper compares the 
ASM with ALM by computing both near and far wake of a Nibe A wind turbine, which combines 
wake velocity, turbulent intensity and vortex structure. Results show that ASM has better prediction 


accuracy and verify it's feasibility on numerical simulation of wind turbine wake. 
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近年 来 , 计算 流体 力学 (CFD) 方法 在 风力 机 流 
场 及 其 气动 性 能 研究 中 得 到 越 来 越 多 的 运用 ， 该 方 
法 能 够 准确 的 描述 出 风力 机 及 其 周围 的 复杂 流 场 . 
常规 的 风机 CFD 计算 , 为 了 保证 风力 机 叶片 及 尾 流 
区 域 的 计算 准确 性 , 需要 建立 相对 复杂 的 实物 模型 ， 
先 将 固体 的 实体 模型 用 专业 的 三 维 建 模 软件 进行 建 
模 , 建 模 完成 后 再 对 整个 流 场 区 域 进行 网 格 划 分 . 三 
维 模 型 的 复杂 性 使 得 在 整个 过 程 中 网 格 的 划分 成 为 
了 一 个 难点 ， 而 且 由 于 在 风 轮 下 游 远 尾 流 区 域内 网 
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格 还 要 保持 一 定 的 精度 ， 网 格 的 数量 大 基 增 加 ， 从 
而 必然 会 导致 计算 量 的 增加 . 

研究 人 员 通 过 将 BEM 理论 与 常规 的 CFD 方法 
相 结 合 , 开发 出 了 致 动 模型 方法 , 即 先 用 BEM 理论 
求解 风 轮 叶片 的 气动 力 ， 再 将 气动 力作 为 体积 力 源 
项 添加 到 N-S 方程 中 求解 , 模拟 叶片 与 流 场 的 作用 
力 。 因 为 模型 中 没有 真实 的 叶片 的 固 壁 边界 ， 所 以 
能 够 大 大 减少 模型 复杂 程度 以 及 网 格 数 量 ， 进 而 节 
约 大 量 的 计算 资源 。 目 前 致 动 模型 主要 有 致 动 盘 模 
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型 、 致 动 线 模型 和 致 动 面 模型 . 

Rajagopalan D 将 致 动 盘 概念 与 CFD 计算 相 结 
fr. 他 采用 时 间 平 均 的 体积 力 分 布 形式 以 及 有 限 差 
分 法 对 风力 机 流 场 的 层 流 N-S 方程 进行 了 定常 计算 。 
沈 翔 和 王 同 光 91 引入 致 动 盘 模 型 , 基于 Garrad Has- 
san 风机 尾 流 实验 , 对 风力 机 尾 流 速度 型 的 分 布 形态 
进行 了 研究 . 许昌 和 韩星 星 O 将 致 动 盘 模型 与 拓展 
k-e i DL: Ein. 引入 汪 流 耗 散 源 项 ， 改 善 尾 流 恢 
复 过 快 问 题 ， 计 算 多 种 风速 下 的 风力 机 尾 流 . 

Sørensen JN 向 和 Shen WZ B 将 致 动 盘 模 型 进 
行 了 改进 , 舍弃 了 致 动 盘 模型 的 轴 对 称 假设 , 提出 了 
用 于 非 定常 计算 的 致 动 线 方法 ,该 方法 更 真实 地 模 
拟 了 风 轮 的 转动 情况 . 卞 风 娇 和 王强 等 9 利用 致 动 
线 模型 ， 基 于 openFOAM 平台 对 风力 机 的 涡 结构 ， 
尾 流 速度 亏损 做 了 数值 模拟 研究 . 

Shen 和 Dobrev 等 人 [9 在 致 动 线 模型 的 基础 
上 发 展 出 了 致 动 面 模型 。 由 于 考虑 了 在 弦 长 方向 上 
受 力 的 变化 ,使 用 致 动 面 模型 计算 出 的 结果 更 加 准 
确 ， 致 动 面 模型 对 咽 尖 涡 以 及 近 尾 流 流 场 的 模拟 精 
度 得 到 了 提高 。 

目前 致 动 模型 的 主要 计算 研究 方向 是 风力 机 叶 
片 的 载荷 特性 ， 转 轮 区 域 流 场 以 及 近 尾 流 区 域 的 流 
场 特 性 . 而 风力 机 远 尾 流 区 域 的 计算 对 于 风 场 开发 、 
微观 选 址 有 着 重要 的 参考 意义 ， 本 文 主要 运用 致 动 
面 模型 对 单 台 Nibe A 型 风力 机 的 尾 流 场 进行 数值 
模拟 计算 , 包括 速度 、 注 流 变化 以 及 涡 量 特性 等 , 并 
与 致 动 线 模 型 的 数值 模拟 结果 进行 对 比 , 验证 致 动 面 
模型 方法 运用 在 风 场 选 址 上 的 可 行 性 . 
1 致 动 面 模型 
1.1 体积 力 的 求解 

致 动 面 模型 实质 上 是 将 风力 机 叶片 简化 成 无 厚 
度 的 平面 ， 在 平面 上 施加 不 连续 的 压力 来 模拟 叶片 
对 气流 的 作用 . 叶片 上 各 点 压力 用 二 维 辟 型 气动 数 
据 , 根据 BEM 理论 计算 确定 . 继而 将 分 布 的 体积 
在 流 场 区 域内 用 三 维 N-S 方程 求解 ,不 可 压缩 N-S 
方程 : 

OV 


m 1 doo 
ou d eMe PEN ERST (1) 


式 中 : /为 作用 在 旋转 叶片 上 的 体积 
如 图 1 所 示 ， 相 对 于 叶片 的 空气 流 相 对 速度 可 
以 从 上 述 速 度 三 角形 中 得 出 : 
Va E V2 x (2r = Vo)? (2) 
AP: 2 为 风 轮 转速 ; V, 为 相对 径 向 速度 ; Vo 为 切 
向 速度 ;el 为 合 速度 ; Va 和 旋转 平面 的 夹 角 (入 


流 角 ) 可 用 下 式 求 得 : 


sepa. e 
$ = tan (ze) (3) 


图 1 叶 素 受 力 分 析 


Fig. 1 The velocity and force in blade element 


当地 攻 角 o = 6— y, REP y 为 安装 角 . 叶片 叶 
展 方 向 单位 长 度 上 的 体积 力 由 下 式 确 定 : 


1 
fap = (L, D) = 5 Vac CLer, Cpep) (4) 


式 中 : Cr = Cla, Re) 为 升力 系数 ; CD = Cp(a, Re) 
为 阻力 系数 ， 是 以 攻 角 和 雷诺 数 为 变量 的 函数 ;er 
为 升力 方向 向 量 ; ep 为 阻力 方向 向 量 ; 雷诺 数 Re 
BIZK c 和 来 流速 度 确定 。 
1.2 满 流 模 型 

由 于 风 轮 在 旋转 过 程 中 ， 对 其 前 后 流 场 产 生 扰 
动 ， 加 快 了 汗 流 动能 的 产生 与 耗 散 ,为 了 尾 流 模 拟 
的 准确 度 ， 需 要 添加 消 流 耗 散 率 源 项 005, 平衡 湛 流 
产生 率 与 耗 散 率 . 文献 [10] 中 研究 了 拓展 k-e dii Ui 
方程 在 水 平 轴 风 力 机 尾 流 计算 中 的 运用 。 消 流 耗 散 
率 源 项 为 : 
Pe 
pk 
式 中 : STU 为 满 流 耗 散 率 ， 用 来 表征 注 流 从 大 尺度 
向 小 尺度 过 渡 时 候 的 能 量 传递 速率 ; RP. 为 消 流 动能 
ERX, 其 中 参数 Ca. W 0.37. 源 项 St 添加 在 风 
轮 前 后 0.5 倍 直径 区 域 .为 更 真实 地 模拟 风机 工作 
时 的 状态 ,本文 添加 了 机 舱 源 项 Sa 和 塔 架 源 项 S. 
具体 如 文献 [4] 所 示 . 本 文 数值 模拟 使 用 拓展 k-e dg 
流 方程 DU, 


2 计算 方法 
2.1 风力 机 参数 


本 文采 用 丹麦 Nibe A 型 水 平 轴 风 力 机 进行 模 
拟 计算 UU, 风机 的 轮 载 高 度 为 45 m， 风 轮 直 径 为 


Se tes (5) 
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40 m, 叶片 设计 采用 气动 经 型 NACA 44 系列 , 额定 
风速 为 13 mst, 根据 G. J. Taylor 等 人 在 该 风机 流 
场 内 做 了 实验 数据 采集 工作 , 以 Nibe A 型 风力 机 为 
模型 在 来 流风 速 为 8.5 msl, 转速 为 3.5 rad/s 的 工 
况 下 , 用 致 动 面 方法 对 风机 尾 流 场 进 行 了 数值 模拟 ， 
并 将 结果 与 实验 测量 值 进行 比较 ，。 
2.2 计算 域 

首先 运用 Gambit 构建 计算 域 并 按照 不 同 区 域 
添加 相应 大 小 的 非 结构 性 网 格 . 为 了 结合 实际 风 场 
中 的 条 件 ， 整 个 计算 域 设计 为 一 个 规则 圆柱 体 ， 半 
径 140 m, 总 长 1203 m。 具 体 划 分 如 图 2. 


图 2 计算 域 划分 示意 图 


Fig. 2 Diagram of computational domain 


2.3 体积 力 的 分 布 

因为 致 动 面 模型 中 不 存在 真实 的 叶片 固 壁 边界 ， 
所 以 风机 叶片 和 流 场 之 间 的 相互 作用 完全 通过 体积 
力 源 项 的 分 布 来 表征 . 通过 上 文体 积 力 的 计算 ,得 
到 的 是 沿 叶片 展 向 单位 长 度 的 体积 力 源 项 ， 未 考虑 
到 真实 情况 下 风 轮 旋转 时 弦 长 方向 上 力 的 分 布 特征 。 
致 动 面 模型 相 比 较 于 致 动 线 模 型 ， 最 明显 的 改进 与 
优化 就 是 在 体积 力 的 分 布 方式 上 考虑 了 叶片 弦 长 上 
的 变化 ， 


图 3 体积 力 分 布 方式 


Fig. 3 Distribution of volume force 
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考虑 到 叶片 弦 长 对 叶片 上 体积 力 分 布 的 影响 ， 
本 文 模型 中 体积 力 的 分 布 采 取 分 段 线性 分 布 方 
式 ba， 如 图 3 HR, URAD 1/4 弦 长 位 置 为 分 
界 点 , 旦 保持 1/4 弱 长 位 置 处 翼 型 的 俯仰 力矩 为 零 。 
这 样 的 分 布 方式 更 符合 真实 情况 下 的 叶片 固 壁 边界 
效应 ,能 有 效 的 改进 近 尾 流 区 域 的 计算 准确 度 . 
2.3 致 动 面 的 辨识 

将 体积 力 源 项 准确 的 添加 到 致 动 面 模型 所 定义 
的 无 厚度 平面 上 是 计算 中 的 重要 环节 . 如 图 4, 叶片 
所 在 平面 即 为 致 动 面 所 在 平面 ，O 点 为 叶片 旋转 中 
心 , P 点 是 风 轮 旋转 平面 内 任意 一 点 ， 任 意 选择 一 
个 叶片 ，@ 为 该 叶片 防线 上 一 点 。 

若 OP x OQ. k > 0, 则 PP 点 在 叶片 弦 线 右边 ; 
反之 则 P 点 在 叶片 弦 线 左边 . 

式 中 , 大 为 Z 轴 正 方向 单位 向 量 (与 来 流风 速 
方向 相反 ). 

若 卫 点 在 弦 线 右边 , H |OP|c/4,(c Junt 
K) 则 PP 点 在 该 叶片 上 ; 若 已 点 不 在 该 叶片 上 , 则 
再 依次 与 另外 两 个 叶片 进行 匹配 ,最 终 可 确定 点 P 
在 平面 上 的 位 置 ， 

叶片 旋转 后 ,OG 随 着 时 间 t 和 角速度 w 变化 ， 
在 每 个 时 间 步 长 上 对 致 动 面 网 格 重复 以 上 步 又 进行 
识别 。 


图 4 致 动 面 识别 示意 图 
Fig. 4 Grid identification of ASM 


通过 上 述 致 动 面 辨识 技术 , 用户 可 以 自 定义 函 
数 UDF 准确 定位 旋转 中 的 致 动 面 位 置 , 这 种 方法 很 
大 程度 上 减少 了 模拟 实际 风电 场 建立 致 动 面 模型 的 
工作 量 ,对 实际 工程 有 很 大 的 实用 价值 


表 1 计算 域 尺 寸 及 网 格 数 


Table 1 Size of calculation domain and number of grids 


区 域 编 号 1 2 3 
长 /半径 (m) 3/140 80/140 3/140 
网 格 数 /104 38 10 53 


总 计 /104 


4 5 6 7 
200/140 200/140 600/140 3/140 
133 17 26 2.5 
281 
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2.4 边界 与 求解 

边界 条 件 方面 , 计算 域 入 口 边界 采用 Inlet 固定 
速度 入 口 条 件 , 出 口 边界 采用 Outlet 零 压力 梯度 出 
口 条 件 , 圆柱 壁面 采用 Symmetry 条 件 . 数值 模拟 采 
用 求解 器 的 默认 设置 ， 压 力 - 耦合 采用 SIMPLEC 
算法 ， 控 制 方程 的 离散 选用 Standard FA, 动量 、 
灌 动 能 和 耗 散 率 均 采 用 二 阶 迎 风格 式 ， 各 松弛 因子 
均 采 用 默认 值 ， 
3 计算 结果 分 析 
3.1 体积 力 的 修正 

致 动 模型 计算 叶片 体积 力 通过 BEM 理论 , 计 
算 各 截面 位 置 的 入 流 角 和 攻 角 大 小 ， 然 后 利用 达 代 
法 计算 诱导 因子 ， 再 根据 二 维 辟 型 气动 数据 查 表 获 
得 每 个 小型 升力 系数 和 阻力 系数 最 后 由 式 (4) 得 出 . 
因此 在 数值 模拟 过 程 中 ,叶片 的 三 维特 性 不 能 够 得 
到 充分 表现 , 且 迭 代 后 拟 合 过 程 也 存在 一 定 误差 , 从 
而 影响 计算 结果 准确 度 . 综 上 考虑 , 添加 一 个 体积 
修正 系数 Ct. 那么 体积 力 可 表示 为 : 


fap = Ctfap (6) 


Cr 拟定 1.0, L1 和 1.2， 分 别 进行 模拟 计算 ， 
选取 风 轮 后 2.5D 及 6D 处 截面 风速 ， 进 行 对 比如 
图 5, 

比较 后 可 见 Cf=1.2 时 风速 曲线 与 实验 数据 更 吻 
合 , 因此 选 定 体积 力 修 正 系数 大 小 为 1.2。 

3.2 计算 结果 分 析 

按照 上 文 所 述 条 件 及 设 定 , 在 FLUENT 中 耦合 
相关 模型 , 对 Nibe A 型 风力 机 的 流 场 进行 数值 模拟 
计算 , 计算 时 间 全 长 为 风机 运转 142 s, 来 流通 过 整 
个 计算 域 . 提取 计算 结果 进行 分 析 . 

图 6 是 风 轮 叶片 平面 处 速度 云图 及 涡 量 图 . 


图 5 两 种 体积 力 修正 系数 对 比 


Fig. 5 Comparison of two correction coefficient 


从 图 6 中 可 以 看 出 ， 风 轮 平面 速度 分 为 3 个 
扇形 ,不同 位置 处 速度 大 小 分 布 清晰 ， 叶 片 外 缘 位 
置 速度 较 大 ， 而 叶 尖 处 速度 最 大 。 从 涡 量 图 可 看 到 
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Zi. 


t9U2 4S LAO -100 


图 6 转 轮 平面 速度 云图 和 涡 量 图 
Fig. 6 Contour of velocity at wind turbine rotor and 


Iso-surface vorticity of ASM 


在 轮 载 高 度 取水 平 截面 分 析 风 力 机 尾 流 场 ， 主 
要 分 析 速 度 与 汪 流 的 变化 , 如 图 7、 图 8 所 示 . 图 7 
中 可 以 看 出 ， 风 机 尾 流 扩张 角 明 显 ， 速 度 变 化 规律 
Hm, 3D 距离 后 趋 于 稳定 , 16~17D 后 风速 基本 恢 
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复 至 来 流风 速 . 由 图 8 可 见 , 油 流 在 风 轮 前 LD 到 风 
轮 后 3D 范围 内 变化 较为 强烈 , 风 轮 后 方 5D 距离 后 
趋 于 稳定 , 在 16~17D 后 基本 可 以 忽略 不 计 ， 这 与 
速度 云图 趋势 基本 一 致 。 
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图 7 水 平面 速度 云图 


Fig. 7 Contour of velocity at horizontal plane 
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图 8 ZK P im zs d 


Fig. 8 Contour of turbulent intensity at horizontal plane 


3.3 致 动 面 模型 与 致 动 线 模型 的 对 比分 析 

考虑 到 实际 风 场 中 风机 间距 一 般 在 5D 以 上 , 如 
果 想 要 将 致 动 面 模 型 运用 于 实际 风 场 的 流 场 计算 之 
中 去 ， 那么 只 对 近 尾 流 区 域 进行 计算 是 不 够 的 . 因 
此 本 文 着 重 对 风机 后 2.5~7.5D 距离 的 风速 情况 作 
研究 分 析 。 在 流 场 中 轮 载 高 度 截 取水 平面 ,分 别 导 
出 在 转 轮 后 侧 2.5D、6D 和 7.5D 三 处 水 平 线 上 速度 
进行 分 析 ， 同 时 将 用 致 动 线 方法 14 模拟 所 得 结果 
一 并 列 出 , 与 实验 数据 进行 比较 , 得 到 图 9. 

从 图 9 中 可 以 看 出 ,三 个 位 置 处 致 动 线 方法 和 
致 动 面 方法 数值 模拟 的 结果 整体 接近 ， 与 实验 数据 
相 比 ， 两 侧 风 速 基本 一 致 ， 最 小 风速 都 出 现在 中 心 
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位 置 . 二 ID 范围 内 尾 流速 度 曲 线 坡度 明显 减 小 , 与 
图 7 水平 速度 云图 相互 验证 . 

由 图 9(a) 可 见 , 在 2.5D 处 , 致 动 面 方法 与 致 动 
线 方法 的 模拟 结果 和 实验 数据 都 很 吻合 ， 整 体 上 都 
呈 V 型 ， 致 动 线 方法 中 心 位 置 速 度 略 有 波动 . 但 从 
图 9(b)、(c) 可 以 看 出 , 在 6D 和 7.5D 处 , 致 动 线 方 
法 计算 结果 和 实验 数据 偏差 较 大 ， 轴 线 位 置 的 实验 
最 小 风速 为 0.7U EA, 而 致 动 线 模型 最 小 风速 只 有 
0.5U。 致 动 线 模型 截面 上 风速 曲线 的 趋 热 和 实验 数 
据 有 较 大 出 入 ,依然 保持 V 型 速度 曲线 ,尤其 风 轮 
位 置 (r/D = 土 0.5D) 风速 明显 偏 小 , 误差 超过 2096, 
而 实验 数据 中 速度 变化 曲线 呈 浅 U 型 . 这 说 明 在 致 
动 线 模型 的 数值 模拟 中 ， 远 尾 流 区 域 动能 耗 散 与 恢 
复 和 真实 情况 误差 较 大 ， 而 致 动 面 方法 的 模拟 结果 
就 相对 准确 很 多 ， 曲 线 走 势 基本 一 致 ， 尤 其 是 在 远 
尾 流 区 域 依 然 保 持 了 较 高 的 吻合 度 . 
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本 文通 过 使 用 致 动 面 模型 对 Nibe A 型 风力 机 
尾 流 场 进行 计算 和 分 析 ， 并 与 致 动 线 模 型 以 及 实验 
数据 进行 对 比 , 得 到 以 下 结论 : 

1) 使 用 致 动 面 模型 结合 CFD 软件 对 风力 机 尾 
流 场 进行 计算 ， 人 简便 易 行 ， 计 算 量 小 ， 模 型 通用 性 
较 强 ; 

2) 与 致 动 线 模型 数值 模拟 结果 对 比 ， 可 看 出 致 
动 面 模型 在 数值 模拟 的 精度 上 有 一 定 优势 ， 尤 其 是 
对 远 尾 流 区 域 的 计算 精度 有 很 大 的 提升 ; 

3) 致 动 面 模型 计算 风力 机 流 场 所 需 计 算 资 源 不 
多 ， 且 对 远近 尾 流 区 域 的 模拟 有 较 高 精确 度 ， 因 此 
可 以 考虑 变 工 况 下 风力 机 的 流 场 数值 模拟 或 者 多 台 
风机 的 计算 . 
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图 9 风 轮 后 2.5D、6D 和 7.5 处 水 平 线 上 致 动 模型 计算 风速 与 实验 数据 对 比 
Fig. 9 Comparison of velocity between ALM and ASM at 2.5,6 and 7.5D behind wind turbine rotor 
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